Elliptic beta integrals and solvable 
models of statistical mechanics 



V. P. Spiridonov 



Abstract. The univariate elliptic beta integral was discovered by the author 
in 2000. Recently Bazhanov and Sergeev have interpreted it as a star-triangle 
relation (STR). This important observation is discussed in more detail in con- 
nection to author's previous work on the elliptic modular double and supersym- 
metric dualities. We describe also a new Faddeev-Volkov type solution of STR, 
connections with the star-star relation, and higher-dimensional analogues of 
such relations. In this picture, Seiberg dualities are described by symmetries 
of the elliptic hypergeometric integrals (interpreted as superconformal indices) 
which, in turn, represent STR and Kramers-Wannier type duality transfor- 
mations for elementary partition functions in solvable models of statistical 
mechanics. 
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1. The simplest elliptic hypergeometric integrals 

In the present paper we discuss relations between a new class of special func- 
tions, called elliptic hypergeometric functions, and solvable models of statistical 
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mechanics. We describe the most complicated known integrable systems defined 
on 2d (two-dimensional) lattices representing continuous spin generalizations of the 
well known Ising model and its various extensions. Actually, these novel integrable 
models correspond to some discretized 2c? quantum field theories. Also we indicate 
connections with the Ad supersymmetric field theories, where elliptic hypergeomet- 
ric integrals have found recently the major application. We start from a brief 
technical introduction to the needed results on special functions and discuss the 
physical systems they apply to in the following chapters. 

General theory of elliptic hypergeometric integrals was formulated in [Sl\ IS31 
IS5j . We skip the structural definition of these integrals and refer for the corre- 
sponding details to a reasonably short survey given in [S9| . 

Let us denote 

oo 

{z;q)oo = Y[(l- zq k ), M<1, zeC, 

k=Q 

the standard infinite (/-product and 

00 1 —1 i+l 

ijio 1 _ zpq 

the standard elliptic gamma function. Below we use the conventions 

T(a,b:p,q) := T(a;p,q)r(b;p,q), T{az ±1 ;p,q) := Y{az;p,q)T{az^ 1 ]p,q), 
T{az ±1 y ±1 ;p,q) := T(azy; p, q)T{az~ 1 y;p, q)T{azy^ 1 ;p, q)r(az~ 1 y~ 1 ;p, q). 

One has the symmetry T(z;p, q) — T(z; q,p) and the inversion formula 

rfo, — ;p,q) = 1, or V (a, a' 1 ; p,q) = — — — - — -, 

which follows from the difference equations 

T{qz;p,q) = 6(z;p)T(z;p,q), T(pz;p,q) = 9(z;q)T(z;p,q), 

where 

8(z;p) = (z;p) 00 (pz~ 1 ;p) 00 
is a theta function. The standard odd Jacobi theta function has the form [WW] 

0i (u|r) = -0 U («) = -J2e viT{k+1/2f e 2wKk+1/2)(u+1/2) 

/cGZ 

= ipi/*e-™(p;p) 0O 8(e 2 ™;p), 

where we denoted p — e 2mT . 

The univariate elliptic beta integral [SI forms a cornerstone of a new power- 
ful class of exactly computable integrals. It is described by the following explicit 
formula 

where T is the unit circle with positive orientation, 

_ (p;p)oo(q;q) 
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and six complex parameters tj, j — 1, . . . ,6, satisfy the inequalities \tj\ < 1 and 
the balancing condition 

6 

(1.2) \{h= P q. 

J'=l 

We use the word "integral" in two meanings. When referred to the exactly 
computable cases, like (ll.lj) or the standard Euler beta integral lying on its bottom, 
it means either the function defined by the left-hand side or, more often, the whole 
identity. In other cases it means an integral representation for a function of interest 
or a class of functions with common structure. 

As shown in |S3j , the left-hand side of relation (jTTTJ) serves as the orthogonality 
measure for the most general known family of biorthogonal functions with the 
properties characteristic to classical orthogonal polynomials (Chebyshev, Hermitc, 
Laguerre, Jacobi, . . . , Askey- Wilson polynomials). In the same paper the elliptic 
beta integral has been generalized to the following function 

IlLrfe^ 1 ;?,?) dz 



(1.3) V(h,...,ts;p, q) = k 



T{z ±2 ;p,q) iz 



where \tj\ < 1 and EL—i *i = {pq) 2 ■ This is a natural elliptic analogue of the Gauss 
hypergeometric function since its features generalize most of the special function 
properties of the 2 -Fi -series |S5l IS9| . For tjt k = pq, j ^ k, F-function reduces 
to the elliptic beta integral and, for this reason, it can be called the elliptic beta 
integral of a higher order. 

In S4] , the author has introduced the following universal integral transforma- 
tion for functions analytical in the vicinity of the unit circle T: 

f dz 

(1.4) g(w;t) = K A(t;w,z;p,q)f(z;t)—, 

Jt lz 

where the kernel 

(1.5) A{t;w,z;p,q) := A(t;w,z) = T(tw ±l z ±l ;p,q), \t\ < 1, 

is a particular product of four elliptic gamma functions. In [SWJ . it was shown that 
this integral transformation obeys the key property making it very similar to the 
Fourier transformation. Namely, its inverse is obtained essentially by the reflection 
t t~ l . 

An explicit example of the pair of functions g{w; t) and f(z; t) in (ll.4[) can 
be easily found from the elliptic beta integral. Indeed, let us denote £5 = tw and 
t 6 = tw^ 1 (so that t 2 n_7=i = Pq)- Then, 

(L6) f{z;t) = r( z ^ P ,q) ' 

4 

(1.7) g(w;t) =r{t 2 ;p,q) J[ T{Ut ]]P , q) J| T^w^-p, q), 

l<i<j<i j=l 

where |iw ±:L |, \tj\ < 1. 

Because of the permutational symmetry, any of the original variables tj can be 
associated with the distinguished parameter t. After fixing t\ = sy,t2 = sy^ 1 and 
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<3 = rx, ti = rx , one can rewrite the elliptic beta integral in the form 

dz 



ip(z)A(r; x, z)A(s; y, z)A{t; w, z)- 

iz 

(1.8) = x(r, s, t)A(rs; x, y)A(rt; x, w)A(st; y, w), 
where rst — ±^/pq and 

= tn7^p,7) = ^b;rtoo( 9 ; 9 )oo^ 2 ;p)^- 2 ;<z), 

(1.9) X(r,a,t) = T{r 2 ,s 2 ,t 2 ;p,q). 

A key application of definition (|1.4p consists in the construction of a tree of 
identities for multiple elliptic hypergeometric integrals with many parameters S4J . 
Using one of the corresponding symmetry transformations, the following relation 
has been derived in 



(1.10) 



x; c, d|£; s) = k / R(c, d, a, b; x, w\s)(f)(w; a, 6|£; s)- 



dw 



iuu 



where the "basis vector" <j> has the form 
(1.11) 



w; a, b\£; s) = T(sa^\ S b^\ J?l w ^^;p, q ), 

V ab 



and the "rotation" integral operator kernel is 
R(c,d,a,b;x,w\s) 



x V sc, sd, 



[H. x -i EH H. 

cd ' V cd ' as ' bs\ pq 



w, 



-w 1 ;p,q 



The function is a generalization of the kernel A(t; x, z), since for ab = pq/ s 2 one 
has the reduction 

pq 

(w;a, — |£;s) = A(s; «;,£). 



as* 



Using the A-kernel, relation (jl.101) was rewritten also in |S8j in a more compact 
form 

f duo 
A(a;x,£)A(P;y,0=K / r(a, /?, 7, S; x, y; t, w)A( T , t, £)A(*; w, £)- ^ 



r(a,P,j,6;x,y;t,w) 



1 



r((5 ±2 , w ±2 ;p, q) 



iw 
±i 



where a/3 = j5 and 



±i 



A(a; x, z)A(/3; y, z)A(f;t, z)A(±; z) rfz 



V 1 o J Jt r(z±^;p,g) iz 

Here we use the condensed notation for parameters of the V- function: V(. . . ax 
■ ■ ■) = V(. ■ ■ ax, ax" 1 . . .). 

The function (j> emerges also in the context of the Sklyanin algebra |Sk) (the 
algebra of the Yang-Baxter equation solutions) , 

SaSp ~ SpS a — i(5o5 7 + SjSq), 



(1.12) 



SoS a — S a So — i - — {SpSj + S-fSp), 
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where J a are the structure constants and (a, /3, 7) is any cyclic permutation of 
(1,2,3). Namely, one has to consider the generalized eigenvalue problems A<f> — 
\B<p, where A and B are linear combinations of four generators S a , a = 0, 1,2,3, 
and A is a spectral parameter. The function <p is defined uniquely up to multipli- 
cation by a constant with the help of two such equations using a pair of Sklyanin 
algebras forming an elliptic modular double |S8j . This algebra represents an elliptic 
extension of the Faddeev modular double [F], but there are actually two different 
modular doubles at the elliptic level which obey different sets of involutions. 

Relevance of the Sklyanin algebra in this setting was noticed first by Rains 
[Rlj . For special quantized values of the parameters, the 0- function reduces to 
the intertwining vectors of Takebe [Tj, which were used by Rosengren in [Ros 
for the derivation of a discrete spin version of relation (|1.10j) . In our case both 
Casimir operators of the algebra (|1.12|) . K = J2 a =o $a an d K 2 = J2a=i JaS 2 , take 
continuous values, i.e. we deal with the continuous spin representations related to 
the integral operator form of the Yang-Baxter equation. 

The following scalar product has been introduced in 



(1.13) (X,^) = « 



X(z)ip(z) dz 
t r(z ±2 ;p,q) iz ' 



It has been shown that both the ^-function itself and the ^-vectors form biorthog- 
onal systems of functions with respect to this measure. In particular, one has the 
relation 



K 



m ±2 ; P ,q) i£ 



(1.14) = . . r -, e ±2i ^, q) VT^S(v - v'), 

where v — cos ip, v' = cos tp' , and S(v) is the Dirac delta-function. (There is a 
missprint in formula (3.2) of |S8j which misses the first factor standing on the 
right-hand side of (jl.l4[l .) Positivity of the biorthogonality measure and of the 
^-function corresponds to the unitarity of representations of the elliptic modular 
double. Setting cd — pq/s 2 , we obtain 



(2/c) 



1 A(s- 1 ;e i¥ '',e ix )A(s;e i¥, ,e ix ) dX 
.! r(e ±2i x;p,q) Vl - ~X~ 2 



(1.15) =r( s 2 , s - 2 ,e ±2iv -p,q) y/l - v 2 5(v - v'), 
where X = cosx. 

The tetrahedral symmetry transformation for 1^-function, discovered in 
can be rewritten in the following form: 

Viax+Wy^^w^^z* 1 ) = T(a 2 , (3 2 , 7 2 , S 2 ;p, q)A(ap; x, y)A( 1 S; w, z) 

(1.16) x VypqP' 1 ^ 1 , y/pqa^y* 1 , s/pqS^w^, Vpqj^z* 1 ) 
= A(«7; x, w)A(aS; x, z)A(/3 7 ; y, w)A(f3S; y, z) 

(1.17) x y(/3x ±1 ,ay ±1 ,(5u; ±1 ,7z ±1 ) 

(1.18) = T(a 2 , /3 2 , 7 2 , S 2 ; P , q)A(a(3; x, y)A(a 7 ; x, w)A{ad; x, z)A(p T , V, «0 

x A(/3S; y, z)A{ 1 6- w, z)V (^/pqa- 1 x* 1 , V\ y/pri' 1 ^ 1 , VM^z^), 
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where a^jS — ±pq. The latter two transformations are obtained by repeated ap- 
plication of the first relation in combination with permutation of the parameters. 
The full symmetry group of the ^-function is the Weyl group W(Et) for the ex- 
ceptional root system E-j |R3j . Therefore, there are 72 = dimW (-EV) / Ss relations 
similar to 1)1.160 . 1)1-17)1 . (11.180 . we just picked up three of them by breaking the 
Ss permutational symmetry and gathering the elliptic gamma functions into the 
A-blocks. 

The outstanding physical application of the elliptic beta integral has been dis- 
covered by Dolan and Osborn [DO] . They have shown that the simplest supercon- 
formal (topological) indices of Af = 1 supersymmetric field theories coincide with 
known elliptic hypergeometric integrals. Exact computability or the Weyl group 
symmetry transformations of such integrals describe the Seiberg duality of M = 1 
theories [Sbj . since they prove coincidence of the corresponding superconformal 
indices. 

In this picture, the left-hand side of the univariate elliptic beta integral eval- 
uation formula 1)1.10 describes the superconformal index of the supersymmetric 
quantum chromodynamics with SU{2) gauge group and SU(6) flavor group. This 
theory has one vector superfield (gauge fields) in the adjoint representation of SU(2) 
and a set of chiral superfields (matter fields) in the fundamental representation of 
SU(2) x SU(6). The elementary particles representing these fields describe the 
spectrum of the theory in the high energy limit, where the coupling constant is 
vanishing due to the asymptotic freedom. In the deep infrared region the theory is 
strongly coupled, all colored particles confine, and one has the Wess-Zumino type 
model for mesonic fields lying in the 15-dimensional totally antisymmetric tensor 
representation of SU{6). The superconformal index of the latter theory is de- 
scribed by the right-hand side expression of formula 1)1.10 . This construction gives 
a group-theoretical interpretation of the elliptic beta integral. After renormalizing 
the parameters tk = {pqY^Vk, the balancing condition takes the form Jlfc=i Vk — lj 
which is nothing else than the unitarity condition for the maximal torus variables 
of the group SU(6). This is the simplest example of the Seiberg duality discovered 
in [Sbj . Further detailed investigation of such interrelations and their consequences 
can be found in [SVlj , where many new elliptic beta integrals on root systems have 
been conjectured and many new supersymmetric dualities have been found. 

The elliptic hypergeometric integrals emerge also in the context of the relativis- 
tic Calogero-Sutherland type models )S7j . However, the first non-trivial example 
of the elliptic hypergeometric functions was found from the exactly solvable mod- 
els of statistical mechanics. Namely, in |FTj Frenkel and Turaev have shown that 
the Boltzmann weights (elliptic 6j-symbols) of the RSOS models of Date et al 
[DJKMO] . generalizing Baxter's eight- vertex model jBaxlj . are determined by 
particular values of the terminating 12 Vn elliptic hypergeometric series (in mod- 
ern notations of |S9j ). The same series has been found by Zhcdanov and the 
author [SZj in a completely different setting, as a particular solution of the Lax 
pair equations for a classical discrete integrable system. In [S2, S3], a family 
of meromorphic functions obeying a novel two-index biorthogonality relation has 
been discovered. It was explicitly conjectured in [S2 that these functions deter- 
mine a new family of solutions of the Yang-Baxter equation for discrete spin models. 
Since the V(t\, ... ,t$;p, q) function is an integral generalization of the latter func- 
tions, in )S6j it was conjectured that the ^-function determines a solution of the 
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Yang-Baxter equation. A simple connection of the terminating i2Vii-series and V- 
function with the Yang-Baxter equation for RSOS models was discussed in |KSj . 
Recently, Bazhanov and Sergeev BS] have shown that the elliptic beta integral can 
be rewritten as a star-triangle relation (STR) which yields a new two-dimensional 
solvable model of statistical mechanics. This is a new important application of 
integral (jl.lj) which is described in the next section. In this paper we show that 
the symmetry transformations for the ^-function have similar interpretation as the 
star-star relations. Moreover, we conjecture that all known exact formulas for el- 
liptic hypergeometric integrals describing the Seiberg duality transformations (at 
the level of superconformal indices) |SVlj . in turn, represent STR and Kramers- 
Wannier type duality transformations }KW|, IWj for elementary partition functions 
in solvable models of statistical mechanics }Bax2| . 

2. The elliptic beta integral STR solution and star-star relation 

In [BS , Bazhanov and Sergeev have interpreted the elliptic beta integral eval- 
uation formula as a star-triangle relation which gave a new solution of this relation. 
In order to describe it, let us introduce the parameter r\ related to the bases p and 
q as 

e" 2 " = pq 

and pass to the additive notation 

z = e iM , x e vx , y ->■ e iy , w e iw . 
Introduce also the exponential form of the parameters 



so that the balancing condition r 2 s 2 t 2 = 1 is satisfied automatically. Finally, denote 

(2.1) W{a; x, u) := A(e Q -"; e ix , e iu ). 
Then relation (|1.8|) can be rewritten as 

S{u\p, q)W{rj — a; x, u)W(a + 7; y, u)W(r] — 7; w, u)du 

(2.2) = x(a,i;P,qW(a;y,w)W(T] -a- j;x,w)W(j;x,y), 
where 

(2.3) S(u; P ,q) = (P;P) ^ g;g)o ° ^(e 2i ";p)g( e - 2i "; g ), 

(2.4) x(<*,r,p,q) = r{r 2 ,s 2 ,t 2 -p,q). 

As observed in [BS , equality (|2.2[) is nothing else than the star-triangle relation 
playing an important role for solvable models of statistical mechanics. It is sym- 
bolized by figure 1 given below, where the black vertex of the star-shaped figure on 
the left-hand side means the integration over u- variable with the weight S(u), and 
M^-weights are associated with the edges connecting the black vertex with white 
ones. On the right-hand side one has the product of three VF-weights connecting 
only white vertices. 

Suppose we have a two dimensional square lattice with spin variables a,b,c, . . . 
sitting at vertices. One associates the self-interaction energy S(a) with each spin 
(vertex). For each horizontal bond connecting spins a and b the energy contribution 
is given by the Boltzmann weight Wf g (a, b), and the energy contribution from each 
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Figure 1. The star-triangle relation. 



vertical bond connecting spins b and d is given by the weight W f g (b, d). The 
variables / and g are called rapidities. Then, as described in detail by Baxter in 
|Bax3l IBax4] , the general STR for these quantities have the following functional 
equations form: 

S(d)W fg (d, b)W fh (c, d)W gh {a, d) = R fgh W fg (c, a)W fh (a, b)W gh {c, b), 

d 

(2.5) ^ S(d)W fg (b, d)W fh (d, c)W gh (d, a) = R fgh W fg (a, c)W fh {b, a)W gh {b, c). 

d 

The second equation is satisfied automatically if the Boltzmann weights are sym- 
metric in spin variables 

W fg (a,b) = W fg (b, a) , W fg (a, b) = W fg (b,a). 

Usually the normalization constants factorize, Rf g h — fghTfg/ffh- Then the weights 
satisfy the unitarity relation of the form 

J2S(d)W fg (a,d)W gf (d,b) = r ~^ffS ab 
d DW 

and the reflection equation Wf g (a, b)W g f(a, b) = 1. 

A subclass of solutions of (|2.5[) emerges from the weights depending only on 
differences of the rapidities, 

(2.6) W fg (a,b) = W(f-g;a,b), W fg (a,b) = W(r)-f + g;a,b), 

where the parameter r\ is called the crossing parameter. Then the precise identifica- 
tion of equality (|2.2p with (|2.5p is reached after setting a = f — g, 7 = g — h (so that 
/ — h = a + 7), equating S(d) to S(u;p, q) and Rf g h to xi a i 1'iPi l) functions, and 
fixing appropriately the range of summation (integration) over the variable d = u. 
We call ([2ll|), (|2l3j), $2A§ the elliptic beta integral STR solution. As show n in [BS] . 
it generalizes many known solvable models of statistical mechanics |Bax2] : the Ising 
model, Ashkin- Teller, chiral Potts, Fateev-Zamolodchikov Zjv-model, Kashiwara- 
Miwa and Faddeev-Volkov models. Moreover, as will be shown below, it comprises 
also a new Faddeev-Volkov type integrable system with continuous spins. 

There is direct relation between spin systems on lattices of three types — the 
honeycomb, triangular, and rectangular lattices. Indeed, one can start from the 
honeycomb lattice, as depicted on the left-hand side of figure 2. Applying the star- 
triangle transformation to each black vertex one transforms the whole honeycomb 
lattice to the triangular one |W] . In a similar way, one can apply STR to each white 
vertex and obtain another triangular lattice having only black vertices. This is quite 
evident and does not require further explanations. However, further transformation 
of the triangular lattice to the square one is more tricky. 
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Figure 2. A honeycomb-triangluar lattice transformation in- 
duced by the star-triangle relation. 



Consider the left-hand side of figure 3. Take the horizontal line in the middle 
of the drawn piece of the lattice. Pick up the triangles above and below it which 
intersect only at one point lying on this line (they are shown in bold lines) . Apply to 
them the triangle-star relation replacing them by stars and continue this procedure 
up and down line-by-line of the resulting lattice. As a result, one obtains eventually 
the square lattice. Taking into account the nontrivial x-multiplier in STR, one 
can thus connect partition functions of the square lattice model to the partition 
functions of two other types of models. 




Figure 3. A triangluar- rectangular lattice transformation in- 
duced by the star-triangle relation. 

In [BS , the parameters x and u in (|2.ip were considered as true spin variables. 
However, because of the x — > —x and u — > —u symmetries, the Boltzmann weights 
W and S depend on their trigonometric combinations. Therefore one can count as 
the true spin variables U — cos it, X — cos a;, etc, with their values ranging from -1 
to 1. The change of the variables in the measure is elementary 

JA> + § = f /(cos " )d " - 2 L m vBw 

The Boltzmann weight W(a;x,u) satisfies the reflection symmetry 



W(a; x, u)W(— a; x,u) = 1, 
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following from the reflection equation for the elliptic gamma function. In terms of 
the spin variables X = cos x and Y = cos y the unitarity relation takes the form 

S(u;p, q)W(r] — a; x, u)W{rj + a; y, u) 



Vi -u 2 



T(e 2a ,e~ 2a ;p,q) 



(2.7) =~^T ^VT^X*6(X-Y). 

S{x;p, q) 



This equality has been established by the author in |S8j . Note that positivity of 
the Boltzmann weights S(u;p, q) and W(a; x, u) corresponds to the unitarity of the 
elliptic modular double representations S8 . In particular, they are positive for 
x,u G [0, 2n], real a such that \y/pqe a \ < 1, and 

1) p* =p, q* = q, or 2) p* = q. 

At the level of superconformal indices, relations similar to (|2.7[) describe the Seiberg 
dualities for gauge field theories with equal number of colors and flavors and the 
chiral symmetry breaking [SV2j . 

Relation (|2.2p is not changed if one replaces W and x by 

W(a; x, u) 



W(a; x, u) 



m(a) 



,„ x m(a)mMm(r) — a — 7) , . 

(2.8) X<*,7;P,g = ( \ \ , ' r X(«,7iP,g) 

m(rj — a)m(r] — 7)771(0; + 7) 

for arbitrary normalizing factor m{a). 

The star-triangle relation is one of the three known forms of the Yang-Baxter 
equation. The second, probably the most popular form, is the vertex type relation 
symbolically written in terms of the i?-matrices as 

(2.9) R( 12 >(A)R( 13 >(A + /i)R (23) (/j) = R( 23 )(/j)R( 13 )(A + /j)R (12) (A), 

where A and /1 are spectral parameters. The third type is referred to as the IRF 
(interaction around the face) Yang-Baxter equation. The star-star relation, which 
was discussed in detail in |Bax3j . belongs to the latter type of equations and has 
the form 



S{g)W x (a, g)W 2 (b, g)W 3 (c, g) W 4 (d, g) 

m(b,c)p{a,b) 
m(a, d)p(c, d) 



(2-10) = R ™?; C ]l { ^ Yl S(9)W{( ai g)W&, g)W£(c g)Wi(d, g), 



<j 

where Wj ;(a, b), Wj(a, b), m(b, c),p(a, b) are two-spin Boltzmann weights and S(g) 
is the spin self-interaction weight (it was omitted in formula (1.1) of [Bax3 ). The 
left-hand side can be interpreted as an elementary partition function for a system 
of four spins a, b, c, d sitting in four square vertices connected by edges to the spin 
g sitting in the square center, and the summation is going over the values of the 
central spin, see figure 4 below. The right hand side has a similar interpretation of 
a statistical sum multiplied by the additional Boltzmann weights associated with 
opposite edges of the square (a, b, c, d). Formula (|2.10[) can be thought of as a 
generalized Kramers- Wannier duality transformation [fCWl IWj . 

Relation (|2.10[) should be compared with the V^-function symmetry transfor- 
mations written in the form (|1.16j) . (|1.17j) . and (|1.18[) . Some of them coincide with 
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Figure 4. A star-star relation for the square lattice. Additional 
Boltzmann weights p and m are indicated by edges connecting 
corresponding vertices on the right-hand side. 



(|2.10l) after appropriate identifications of the Boltzmann weights. For instance, 
equation (| 1 . 1 6[) corresponds to the choice 

W 1 {a,g) = A(a;x,g), W[(a, g) = A^/T 1 ; x, g), 

W 2 (b,g) = A(fi;y,g), W^(b,g) = A(- s /pqa~ 1 ;y,g), 

W 3 (c, g) = A( w, g), W 3 '(c, g) = A(<T 1 ; w, g), 

W 4 (d,g) = A(^pq6;z,g), Wi(d,g) = A( 7 -\z,g), 

where af3^8 — 1, g is the integration variable for the V- function, and S(g) — 
n/V(g ±2 ;p, q). Other factors have the form 

r(a 2 ,/3 2 ;p,g) 
Tij-^S- 2 ^^) ' m(b,c) = m(a,d) = l, 

(2.11) p(a, b) = A(aj3; x, y), p(c, d) = A{a/3; w, z). 

A similar interpretation is valid for relation (j!.17|) . It corresponds to the choice 

W 1 (a,g) = A(a;x,g), W[(a,g) = A(/3;x,g), 

W 2 (b,g) = A(^pq T ,w,9), W^{b,g) = A(^pq5;w, g), 

W 3 (c,g) = A(/3;y,g), W^{c,g) = A(a;y,g), 

W 4 (d,g) = A(Vpq<5;z,g), Wi(d,g) = A{y/pqr,z,g), 

where, again, af3^5 = 1 and S(g) = K/T(g ±2 ;p,q). As to other factors, R = 1 and 

m(b, c) = A(^/pq(3-y; y, w), m(a, d) = A(^/pqf3j; x, z), 

(2.12) p(a,b) = A(^/pqa^;x,w), p(c, d) = A(^/pqaj; y, z). 

There are three star-star relations for the Ising type models listed in [Bax3 as 
equations (2.16), (5.1), and (5.2). Our first option (|2.11[) corresponds to relation 
(5.2) in [Bax3 . Relations (2.16) and (5.2) in |Bax3j are obtained from each other 
by a reflection with respect to the lattice square diagonal (b, d). Our second option 
(|2.12[) corresponds to relation (5.1) in |Bax3j with nonconstant p- and m-weights. 
However, we have the third nontrivial form of the symmetry transformation for 
the V- function (|1.18j) . It corresponds to a more complicated type of the star-star 
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relation 

(2.13) S (9)Wi (a, g)W 2 (b, g)W 3 (c, g)Wi(d, g) 



R 



m(6, c)p(a, b)t(a, c) 



E 

g 



S(g)W{(a,g)Wl(b,g)W£(c,g)Wi(d,g), 



m(a, d)p(c, d)t(b, d) 
where t(a, c) is a new diagonal Boltzmann weight. Explicitly, we have 



Wi(a,g) = A(a;x,g), 
W 2 (b ) g) = A(/3;y ) g), 
W 3 {c,g) = A(- s /pq'y;w,g), 
W 4 {d,g) = A(^/pqS;z,g), 



W{(a, g) = A(y/pqa 1 ;x,g), 
Wl(b,g)=A(^pqp- 1 ;y,g), 
Wj t {c,g) = A( 7 - 1 ; M , g), 
Wi(d,g) = A(S- 1 ;z,g), 



where a/?7<5 
and 



1. Other factors in (|2~T3)l arc R = r(a 2 , /3 2 ;p, q)/T{-y- 



;p, q) 



m(b, c) 
p(a,b) = 
t(a, c) = 



= A(^/pg/37;y,w) 
: A(af3;x,y), 
A(^/pqa-f;x,w), 



m(a, d) = A(y/pq/3j; x, z), 
p(c, d) = A(a/3; w, z), 
t(b,d) = A(y/pqaT,y,z). 



Perhaps, this type of the star-star relation was not considered in the literature 
before. Note that all such relations represent symmetry groups of the partition 
functions. In the case of ^-function this is W(E?), i.e. one has much bigger 
symmetry than that seen explicitly in the chosen spin system interpretation. We 
have described thus a new (elliptic hypergeometric) class of solutions of the star-star 
relation which should lead to new solvable models of statistical mechanics similar 
to the checkerboard Ising model. Known systems of such type were investigated 
in detail in BSt] . A natural general conclusion from our consideration is that the 
symmetry of STR can be richer than a direct sum of symmetries of the Boltzmann 
weights and the lattice. 

3. A hyperbolic beta integral STR solution 

We describe now another solution of the star-triangle relation associated with 
the modified form of the elliptic beta integral when one of the bases p or q can lie 
on the unit circle [DS2]. It simplifies also consideration of the degeneration limits 
to q-beta integrals of the Mellin-Barnes type (hyperbolic beta integrals). 

First we describe the modified elliptic gamma function introduced in S3] . It is 
convenient to use additive notation and introduce three pairwise incommensurate 
quasiperiods u)±, ll>2, <^3 together with the definitions 



2m- 



p = e 
p = e 



2m- 



r 

f ■ 



2m- 



2m^ 

a 3 



(3.1) q = e 

Here q, p, and f are particular (r — > — 1/t) modular transformations of q,p, and 
r. Assume that Im^Wi/wa), hxi^s/cJi), Im^a/u^) > 0, or \q\, \p\, |r| < 1. Then the 
modified elliptic gamma function is constructed as a product of two elliptic gamma 
functions 



T(e 27rl "2 ;p, q)T(re 

e -Tr£ 3 ,3(«;w)p(e-27ri 



(3.2) 



2ni 



\r, q) 
3 ;f,p), 
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where -63,3(11; a;) is the third diagonal Bernoulli polynomial (for the general defini- 
tion of such polynomials, see Appendix A), 



B. 



3,3 




2 / W1W2W3 

n— 1 ' 



The G(u; w)-function satisfies the following system of three linear difference equa- 
tions of the first order 

G(u +wj;w) = ^e 27 " 1 ^ ' ;p)G(u;u), 
G(u + uj 2 ; lj) = 9(e 27rl ^ ; r)G{u; u), 
G(u + lj 3 ; u) = e-* {B ^ u ^G{u- u), 
where B 2 ^ 2 (u;lu) is the second diagonal Bernoulli polynomial, 

U 2 U U U)t U) 2 1 

B 2 , 2 (u; w) = + + 7^+o- 

The second equality in Q3.2p follows from the fact that both expressions for G(u; u) 
satisfy the above set of equations and the normalization G(^ Y^k=i w fc;<^) = 1< 

It is easy to see that G(u;uj) is well defined for |p|,|r| < 1 and |g| < 1, the 
\q\ = 1 case being permitted in difference from the T(z;p, g)-function. Evidently, 
we have the symmetry relation 

G(u; u)x,u}2,u)a) = G(u\ u> 2 , wi, w 3 ) 

and the reflection equation 

3 

G(a;u)G(b;u>) = 1, a + 6 = ^w fe . 

k=l 

For Im(aJi/u; 2 ) > 0, we can take the limit lu 3 — > 00 in such a way that 
Im(w 3 /wi), lm(oj 3 /u} 2 ) +00 

and p, v — > 0. Then, 

(3.3) lim G(u; uj) = j(u: 1^2) = — ' , °° ■ 

For Re(wi), Re(w 2 ) > and < Re(it) < Re(wi + uj 2 ) this 7-function has the 
following integral representation 

e ux dx" 



(3.4) 7(u;wi,w 2 ) = exp 1 , . q (1 _ e u, ix)(1 

which shows that j(u; Wi,u; 2 ) is a meromorphic function of u even for uji/lj 2 > 0, 
when |g| = 1 and the infinite product representation Q3.3[) is not valid any more. 
The inversion relation for this function has the form 

7(«;w 1 ,w a )7(wi+W2-«;wx»W2) = e 7 ^ 2 ' 2 ^. 
For more details on this function see Appendix A. 

Let Im(wi/w2) > and Im(w3/a;i), Im^a/o^) > 0, and let six complex param- 
eters gfc, k = 1, ... ,6, satisfy the constraints Im(gk/ui 3 ) < and 

6 

(3.5) y^fffc = LOl + LQ 2 + UJ 3 . 

k=l 
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Then [DS2] . 

J ,,./ 9 Cr ±2U W - LJ - 

17 -"a/^ v. ' ' l<fc<;<6 

where the integration goes along the straight line segment connecting points — W3/2 
and W3/2, and 

(3.7) «= - 2 ^(g;g)^ 

(q; q) 00 (p;p) 00 (r-^r)^ ' 

Here and below we use the shorthand notation 

G(a ± 6; w) := G(a + b, a - 6; w) := G(a + b; uj)G{a - b; u). 



The proof of equality (|3.6|) is rather simple. It is necessary to substitute in it 
the second form of G(u; w)-function (|3.2|l . check that all exponential factors cancel 
and, after a change of notation, the formula reduces to the standard elliptic beta 
integral. 

Let us introduce the crossing parameter 



3 



1 3 

fc=i 

and denote 

(3-8) g h2 = -a±x, g 3A = a + 7 - <rj ± y, g 5fi = --f±w, 



so that the balancing condition (|3.5p is satisfied automatically. Introduce also the 
modified Boltzmann weight, or the kernel for the modified form of the integral 
transformation (11.41). 



W (a; x, u) = G(a — 77 ± x ± u; ui). 
Then relation (|3.6|) can be rewritten as 

r w 3 /2 

4>(u; u)W'(r] — a; x, u)W'(a + 7; y, u)W (ry — 7; iu, w)ohi 

-W3/2 

(3.9) = x(", 7; u)W'(a; y, w)W'(j] - a - 7; x, w)W'(j; x, y), 



where 



' ig— TiB a ,2(2u;un,ii;3)g/ e — 4iriit/«a . p)0( e -4 '"' in / LJ i • r ) 



' kG(±2u;u;) k 

(3.10) 7; w) = G(-2a, -2 7 , 2a + 2 7 - 2?;; w). 

Substituting the second form of the modified elliptic gamma function, we find 

rrr// \ / 4?ri / . , 3a(a; 2 + m 2 ) \ \ 
W(a;aj,u) = exp — — B3.3 (a - 77; w + — - '-) 

(3.11) xA (e 2 ^-^,-;pA 

We see that this Boltzmann weight is obtained from (12.1[) after a reparametrization 
of variables and multiplication by an exponential of a quadratic polynomial of the 
spin variables. This means that there exists a nontrivial symmetry transformation 
of the star-triangle relation modifying its solutions in the described way. 
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The distinguished property of the modified elliptic beta integral is that it is 
well defined for \q\ = 1. Therefore the limit 0J3 — > 00 leads to g-beta integrals well 
defined in this regime as well. Let Re(wi), Ref^) > 0. Then, for 0J3 —> +ioo, 
one has p, r — > and G{u\u>) goes to 7(1*; u}±, W2)-function. Let us substitute 
g$ = ^. j LJk — A in formula (I3.6[) . where A = Y], g^, and apply the inversion 
formula to the corresponding modified elliptic gamma function. Then the formal 
limit W3 — > +ioo reduces this integration formula to 

Il5=i7(0*±u;w) . _ (g; g)oo Ilt<i<*<B 7Cffj +dk;u) 

— —au = —ZLU2-, s =r — = ! 

i{±2u,A±u;uj) (g;g)oo IT k =ii( A - 9k;u) 

where the integration contour is the straight line for Re(<7fc) > or the Mellin-Barnes 
type contour, if these restrictions for parameters are violated. Let us remind also 
that 

(q; q)oo V w 2 

where V~— I = e -7 "/ 4 since for = ia, a > 0, the square root should be positive. 

Let us introduce parameter g$ anew (it should not be confused with the previous 
variable ge which we have eliminated) using the condition 

6 

(3.13) fffc = ui + L02 

fe=l 

(note the difference with (|3.5[1 ). Now we can apply the inversion formula to 7- 
functions to move some of them from the denominator of the integral kernel to its 
numerator. It is convenient here to define the hyperbolic gamma function 7^) (it): 

(3.14) 1 (2) {u;uj) = e-^ B2 ^ u ^j(u;uj). 

Then, after the change of the integration variable u = iz, the integral p,12[) takes 
the compact form 

J -°° 7 \^ 1Z ^> l<j<k<6 

Validity of the described limit W3 -> 00 at the level of integrals was rigorously 
justified in [R2J using a slightly different notation. Integral (13. 15)) was proven first 
(using a different approach) by Stokman |Stj who called it the hyperbolic beta 
integral. We followed the formal limiting procedure suggested in |DS2] . 

Similar to (|3.8|) . let us fix the parameters as 

(3.16) g 1<2 = -ol ± ix, 33,4 = a + 7 - i] ± iy, g 5fi = -j±iw 

with the crossing parameter 77 = — (wi+W2)/2. Then formula (|3.15j) can be rewritten 
as a star-triangle relation 

/oo 
S{z)W{rj — a; x, z)W(a + 7; y, z)W{rj — 7; w, z)dz 
-00 

(3.17) = x(a, "f)W(a; y, w)W(t] - a - 7; x, w)W(r, x, y), 
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where 

W(a; x, z) — 7^ 2 -* (a — i] ± ix ± iz; ui), 



i 2 sinh 2*2 sinh ^ 

S(z) - 



2y/u> 1 u! 2 J (2 H±2iz; uj) ^^\ui 2 
(3.18) x(a,7) =7 (2) (-2a, -2 7 , 2a + 27-27;; cj). 

These Boltzmann weights are positive for real x,z,a, r\ < a < —77, 77 < 0, and 
either real or u>l = ui 2 . 

Consider a particular reduction of integration formula (I3.15[) . For this we re- 
place parameters gj — > gj + i/i, j = 1, 2, 3, and gj — » gj — i/j, j = 4, 5, 6. Since the 
integrand is symmetric in z we can rewrite the left-hand side as 



00 lfj=i 7 (2) (gj + i/J ± iz, ffj+3 - i/x ± iz; uj) 

7 ( 2 )(±2iz;w) Z 

II ^ (2) ~ iz ' &'+ 3 + iz ; u) Pl (z)p 2 (z)dz, 







where 



e -2 7 r(x+ Al )(w 1 1 +w 2 x ) 

Pi (2) 



and 



7 ( 2 )(±2i(z + ui) n-^+oo 
3 

p 2 ( z ) = e M^)(^+^) "Q 7 (2) (5j + 2iAi + iz? ffi+3 _ 2iM _ iz . w) 

j=i 

— S— I -2^(w!+w 2 )+^ ]T?=i ( 9j 2 + 3 -Sj+(9i-9j + 3)("l+"2) 

— > e \ v / / (1 + o(l)J. 

On the right-hand side we find 

3 6 
3=1 fc=4 

where 

^3(3) = ]J ^ (2) + 5 fe + 2i ^' -9j+3 + .9fc+3 - 2i/j; w) 

l<j<fc<3 

¥23i<,-< fc <3 ( B2,2(9j+3+9fc + 3-2i/i)-B2,2(Si+9fc+2i/i) I 

-> e - - V / (l + o(l)). 

One can check that the leading asymptotics of p$(g) coincides with that of the 
P2 (^-function. Taking the limit (i — » +00, which is uniform, one comes to the 
following exact integration formula [B] 

3 6 



J[ 7 (2) {g 3 ~ iz, g j+3 + iz; uo)dz = ^UT 2 J] JJ 7 (2) (<?j + g k ;u), 

-°° 3 = 1 j=l k=A 

where YX=i 9k = ui + u 2 . 
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Let us change notation for the integral parameters 

9i = -a + iar, g 2 =a + ^-r}-\- iy, g 3 = -7 + iw, 
94 = -a-ix, g 5 = a + 7 - 77 - iy, g% = -7 - iw, 

where r\ = — (w!+cj 2 )/2 is a crossing parameter; the balancing condition ^fe=i 9k = 
—2r/ is satisfied then automatically. Now we can rewrite equality (|3.19l) as the star- 
triangle relation 

f°° dz 
/ W{rj — a; x, z)W(a + 7; y, z)W(r) — 7; w, z) 

(3.20) = x W(a; y, w)W(t] - a - 7; x, w)W(-f; x, y), 
where the Boltzmann weight is defined as 

W(a; x, z) — [a — 77 ± i(x — z);u>) 
and the normalization constant is 

(3.21) x = 7 (2) (-2a, -27, 2a + 2 7 - 2??; u>). 

Note that W(a;x, y) = W(a;y,x) and W(a;x, y) > in the same domain of 
parameters as before. Denoting cj = b, 0J2 = b^ 1 , and a = —(6 + b~ 1 )8 / (2ir) one 
can see that W(a; x, z) coincides with the Boltzmann weight of the Faddeev-Volkov 
model [FVl |VF] de noted a s We (x - z) in |BMS] (our differs by sign from the 
definition chosen in |BMS|) up to some normalization factor Fg. 

We thus see that the Faddeev-Volkov model solution of the star-triangle relation 
|VF] is a particular case of our hyperbolic beta integral STR solution (|3 . 1 8|) . The 
fact that the left-hand side of STR for the Faddeev-Volkov model represents a 
particular limiting case of the elliptic beta integral was known to the author already 
in 2008. After seeing BMS and understanding this fact, the author was interested 
whether a similar interpretation exists for the elliptic beta integral itself. However, 
this idea was not developed further, partially because the origin of the normalizing 
factor Fe given in [BMS] was not understood at that time. Fortunately, Bazhanov 
and Sergeev have independently answered this question in [BS . 



4. Partition functions 

The partition function of a homogeneous two dimensional discrete spin system 
on the square lattice with the Boltzmann weights W(a\ m, Uj) (|2.ip and S(uj) (|2.3p 
has the form 

— a; tife, ui) J^J S(u m )du m , 

where the first product is taken over the horizontal edges (ij), the second product 
goes over all vertical edges (k,l), and the third product (in to) is taken over all 
internal vertices of the lattice. Let us take the elliptic beta integral STR solution 
of [BS] and consider the contribution to Z coming from a particular vertex u 
surrounded by the vertices u\, u 2l U3, 114: 

I S(u)W(a; Ui, u)W(a; u, Uz)W(r) — a; 112, u)W(ri — a; u, u±)du. 
Jo 
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Substituting explicit expressions for the weights, one can easily see that this integral 
is equal to the elliptic hypergeometric function V(ti, . . . ,tg;p, q) described above 
(jl.3[) with the following restricted set of parameters 

*2, *3, U} = { e M 'e 27dUl , e^e" 2 ™ 1 , e tt -"e 2witl3 , e Q -"e~ 2Tritl3 }, 

{t 5 ,t 6 ,t 7 ,t 8 } = {e~ a e 2viU2 ,e- a e~ 27riU2 ,e~ a e 2niUi ,e- a e- 27riUi }. 

In total, there are 5 independent parameters, instead of 7 for generic ^-function 
(in addition to the bases p and q). Therefore we conclude that the full partition 
function Z is given by an elliptic hypergeometric integral constructed as a tower 
of intertwined (restricted) elliptic analogues of the Gauss hypergeometric function 
similar to the Bailey tree for integrals |S4j . 

According to the general reflection method used in [BSj . the leading asymp- 
totics of the partition function for two-dimensional N x M lattice when its size goes 
to infinity, N, M — > oo, has the form 



Z = m(a) 

Af,M-s-oo 



A M 



where m(a) is the normalizing factor for Boltzmann weights which guarantees that 
on the right-hand side of STR the \- multiplier (|2.8[) is equal to unity, X = 1. This 
condition is satisfied if 

(4.1) , r(e- 2 ";p,g) = l, or m(a + rj) = r(e 2a ;p, q)m(-a). 

m(r] — a) 

Let us introduce the function 
(4.2) 

M(x;p,q,t)=e X p( V — ( ^ } " ^ _ T(xt^/pqt;p,q,t 2 ) 

where 



, (i , ..(l-p n )(l-q n )(l+t n )J V^^p-ql^^^ 2 )' 



T(z;p,q,t) = J] (l-^W)(l-^V+ l p*+V +1 ), H,kl,|i| < 1, 
j,k,l=0 

is the second order elliptic gamma function satisfying the t-difference equation 

T(tz;p,q,t) = T(z;p,q)r(z;p,q,t). 
The reflection equation T(z~ 1 ;p, q, t) — T(pqtz;p, q, t) leads to the equality 

M(x' 1 ;p,q,t)M{x;p,q,t) = 1. 
It is easy also to check validity of the functional equation 

M(x- 1 p,q,t)M(t~ 1 x;p,q,t) = T^x^^^p, q^j , 

which is equivalent to (|4.ip after fixing t = pq and x — e 2a . Therefore we find the 
needed normalizing function 

(4.3) m(a) = M(e 2a ;p,q,pq), m(a)m(—a) = l. 

The function — log m{a) defines thus the free energy per edge of the integrable 
lattice model under consideration. 
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Now we discuss the partition function for the general hyperbolic beta integral 
solution of the star-triangle relation (|3.18j) . The needed normalization constant 
m(a) is found from the equation 

(4.4) m( - a ) x 7 ( 2 )(-2a;a;) = l, or m(a/2 + 77) = 7 (2) («; w)m(-o/2), 

771(77 — a) 

where 77 = — {lji + w 2 )/2. 

Let us define the function 

7 (3) (w+ I Lfc=i w fc + w 3 ;wi,w 2 ,2u; 3 ) 
M(«;wi,w 2 ,a;3) = 3 S - ^ - ' 

where 7^ 3 ^-function is the hyperbolic gamma function of the third order defined in 
Appendix B. Using the integral representation for it, we can write 

e vx dx^ 



( ' Kxa f 
(4.5) n(u; u)) = exp I / 

v JR- 



i0 (e" i:c - l)(e^ x - l){e^ x + 1) x 



where v = u + Yl\=i w fc/2 a-nd 



= £3,3(1" + w 3 ;wi,w 2 ,2w 3 ) - £3,3(1;; wi,cd2,2a;3) 



3^ 3 2 



2wiW2 V 12 
For a special choice of the third quasiperiod variable w 3 = uj\ + uj 2 , this function 
appeared for the first time in |LZj . 
Using the reflection equation 

3 

7 (3) (^w fc - u; wi,w 2 ,a; 3 ) = 7 (3) (w; u>i, u> 2 , w 3 ) 

and the difference equation 

7 (3) (it + w 3 ; £Ji , w 2 , w 3 ) = 7 (2) (it; Wi, w 2 ) 7 (3) C"; 0J2, W3), 
one can easily check that /i(u; w)/i(— ii; u) — 1 and 

fj,(u;u)n(u - w 3 ;u;) = 7^(71 + -{u x + lo 2 - w 3 ); Wi, w 2 ). 

The latter relation coincides with equation (4.4) for u = 2a and w 3 = U)\ + uj 2 . 
Therefore we find the free energy per edge as — logm(a), where 

(4.6) m(a) = n(2a; u>i, u> 2 , w± + uj 2 ). 

By construction this function satisfies also the reflection equation 77i(a)77i(— a) = 1. 
Denoting u>i = b,uj 2 = b^ 1 and substituting the infinite product representation of 
the 7^ 3 ^-function given in Appendix B, we find the expression 

771(a) = exp ^-Tria 2 - |^(1 - 2(6 + 6" 1 ) 2 ) 

(ge 2 ™/ b ;q 2 )oc fr 1 + e ""/( fc+fc ~ V +1 g 2fc 

( X ( ge 27riub. ^2)^ l-[ o l- e Triu/(b+b-i)pj+lq2k> 

where it is assumed that \q\ < 1, q = e 27ni>2 , q = e ~ 2m / b2 , an d p = e - 7I " 1 /( 1 + fc2 ). 

We turn now to the Faddeev-Volkov model solution of STR (|3.20p . In this case 
we have no self-interaction of the spins sitting in lattice vertices, and the Boltzmann 
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weights attached to edges are simplified. But the partition function asymptotics 
is the same as in the previous case, since evidently the normalizing constant 771(a) 
is found from the same equation (14.41) . The free energy per edge for this model 
was computed already by Bazhanov, Mangazeev, and Sergeev in BMS], where the 
Boltzmann weights normalizing factor was denoted as Fg . Comparing this constant 
with our 777(a), we see that they coincide for a — —(b + b~ 1 )8/(2ir), as necessary. 
However, our infinite product representation of 777(a) in (|4.7|1 differs drastically from 
that given in [BMSJ (which was the source of author's old time confusion). 



5. Conclusion 

After the discovery of elliptic hypergeometric integrals, for a long time the au- 
thor was drawing attention of experts (including the second author of |BSj ) in two- 
dimensional conformal field theory and solvable models of statistical mechanics for 
a potential emergence of such functions in these fields. The connection between the 
elliptic beta integral and the star-triangle relation found in jBSj and the star-star 
relation described above confirms this expectation. However, the nature appeared 
to be much richer than it was imagined in |S2, S6, S9J. As mentioned already, the 
Dolan-Osborn discovery of a stunningly unexpected coincidence of elliptic hyper- 
geometric integrals with superconformal (topological) indices in four dimensional 
supersymmetric gauge theories strongly pushed forward the development of the 
theory and raised many interesting open questions [DOl [SVl] . The interpretation 
of exact computability of the elliptic beta integrals as the confinement phenome- 
non in quantum field theory is a new type of conceptual perception of the exact 
mathematical formulas. 

As to the models considered in this paper, we have described a generalization 
of the Faddeev-Volkov solution of STR [VFj with the continuous spin variables 
taking values on the real line, which was not considered in BS . It has some 
nontrivial self-interaction energy for each vertex and a more complicated form of 
the Boltzmann weights for edges, though the free energy per edge appears to be 
the same as in the Faddeev-Volkov model. In }VF] the Yang-Baxter equation was 
proved using the quantum pentagonal relation |FK V] . It would be interesting to 
interpret in a similar way the model we have described here. Some time ago the 
author has tried to find an elliptic analogue of the pentagon relation in analogy 
with the constructions described in [V], but could not do it yet. Clearly the elliptic 
beta integral gives already an analytic form of that wanted operator relation, but 
it is hard to formulate it in terms of the commutation relations of some explicit 
operators. 

In |FVj . Faddeev and Volkov have considered a lattice Virasoro algebra and 
described an integrable model in the discrete 2d space-time (it was discussed also 
in detail in |FK V] ) . The elliptic beta integral yields more general solutions of 
STR than that of jVFj . and it is natural to ask for explicit realization of the 
corresponding models similar to |FV| . During the work on |SV1] . G. Vartanov 
and the author have suggested that there should exist some elliptic deformation of 
the primary fields V a (z) built from free 2d bosonic fields (in the spirit similar to the 
situation discussed in [SWy| ) such that the three point correlation function would 
be given by the elliptic beta integral and the four point function would be described 
by the lAfunction satisfying the elliptic hypergeometric equation |S9j (so that the 
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tetrahedral symmetries of the V- function would describe the s-t channels duality). 
Unfortunately, such hypothetical results are not conceivable at the present moment. 

From the point of view of superconformal indices the partition function as- 
sociated with the elliptic beta integral solution of STR looks as a superconfor- 
mal index for a particular SU (2)-quiver gauge theory on a two dimensional lat- 
tice. Recently there was a great deal of activity on interrelations between four- 
dimensional super- Yang-Mills theories and two-dimensional field theories, see, e.g., 
[AGTl ICNVl IGPRRl INSl |SWy| . In this framework, the elliptic hypergeomet- 
ric integrals describing superconformal indices of Af = 2 quiver gauge theories have 
been interpreted by Gadde et al in [GPRRJ as correlation functions of some 2d 
topological quantum field theories. 

Therefore it is natural to expect that superconformal indices of all four di- 
mcnsionanl CFTs are related to discretizations of 2d CFT models and other inte- 
grable systems. A connection of the Yang-Baxter moves with the Seiberg duality 
has been briefly discussed in |HV| . In this context, superconformal indices of all 
quiver gauge theories should correspond to full partition functions of some spin sys- 
tems. In view of the abundance of supersymmetric dualities and rich structure of 
the corresponding superconformal indices (twisted partition functions) [SVT], the 
author considers the present moment only as a beginning of uncovering new two- 
dimensional and higher-dimensional integrable models hidden behind the elliptic 
hypergeometric functions. 

For instance, the elliptic Seiberg integral defined on the BC n root system reads 



[D5I1 159) : 



(5.1) 




where \t\, \t m \ < 1, t 2n 2 J]l = i t m = pq, and 

_ (p;p)£,(g;g)£, 
(47r)"n! 

After some work, this formula can be given the STR type shape 




S(u; t, p, q)W(r) - a; x, u)W(a + 7; y, u)W{r] - 7; w, u) [du], 



(5.2) 



W t (a; y,w)Wt{rj -a- 7; w, x)W t {j- 1 x, y) 



where we denoted 



n 



T(t;p, q)dv,j 
r(**;p, g) 



[du] = K n Y[ - 

and the crossing parameter rj is defined as 

e- 2 " = p 



pqt 



The Boltzmann weights have the form 



T(te ±iu ' ±iah ;p,q) 
r(e ±kl J ±iu <=;p, q) 



n 



1 



(5.3) 



S{u;t,p,q)= J] 



n 



T(e ±2 ' m i;p,q) 



l<j<k<n 



J'=l 
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and 

1 - 

(5.4) W(a;x,u) := T\r(^e a e ±ix e ±i ^; P ,q), 

to a) J -- L 

(5.5) W t (a;x,y) := — —T\r( v /pqe a t 3 ~ I1 i 1 e ±ix e ±iy ;p,q) 1 

TOla) J - J - 

and satisfy the reflection relations 

W{a; x, u)W(-a; x, u) = 1, W t (a; x, y)W t (-a; x, y) = 1. 

The normalization constant m{a) for n > 1 has a substantially more complicated 
form than that for n = 1. To describe it we introduce the function 

r(xts 2 ,xt 1 -"s;p,g,t,s 2 ) 

M{x;p,q,t,s) = — — — ^ 

^(xtSjXt 1 n s' 1 ;p,q,t,s z ) 

a particular ratio of four elliptic gamma functions of the third order. More precisely, 
one has 



T(z;p,q,t,s) := JJ 

i,j,k,l=0 



1 _ z -l pl+ l qJ + l t k+l s l+l 
1 — zp l qH k s l 



for zeC*, |s| < 1, with the reflection equation T(z,pqtsz 1 ;p, q, t, s) = 1 

and the difference equation Y(sz;p, q, t, s) = T(z;p, q, t)T(z;p, q, t, s). Then, 

(5.6) m(a) = M{e 2a ;p,q,t,pqt n - 1 ), 

with the standard reflection relation m(a)m(—a) = 1. 

Let us discuss a physical meaning of the obtained model. Consider a honeycomb 
lattice on the plane with two types of vertices - black and white with two adjacent 
vertices always being of different color, see the left-hand side of figure 2. Into 
each white vertex we put an independent single component continuous spin x. 
Into each black vertex we put n independent spins Uj, j = 1, ...,n, or one n- 
dimensional spin with n continuous components. These "spins" are quite different 
from those of the Ising model where they take only the values +1 and —1 (i.e., they 
represent the fields and not the compact spins). In different words, one associates 
to each black vertex the S'P(2n)-group space related to the root system BC n . We 
associate with each black vertex the self-interaction Boltzmann weight (|5.3p with 
an additional interaction between "spin" components in the internal space. To each 
bond connecting "black-and-white" vertices we attach the Boltzmann weight (|5.4|) . 
Then, on the left-hand side of (|5.2j) we have the partition function of an elementary 
cell with the black vertex in the center and the integral taken over the w^-spin 
values. If we apply this star-triangle relation to each black vertex we come to a 
different spin system associated with the plain triangular lattice having only the 
white vertices with the bond Boltzmann weights described by the function (|5.5p . 
see the right-hand side of figure 2. Such a transformation of lattices looks quite 
similar to a transformation of the honeycomb-triangular Ising systems considered in 
[Wj . Perhaps there exists also another STR type duality transformation involving 
only the white vertices (with some self-interaction) allowing for a transition to yet 
another triangular lattice system. In addition to this uncertainty, it remains also 
unclear the free energy per edge of which model is described by the function (|5.6[) . 

Positivity of the Boltzmann weights of this model can be analyzed along the 
lines of elliptic modular double involutions discussed in |S8j . In particular, these 
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weights are clearly positive for x,Uj £ [0,27r], real t,a and p := \y/pqe a \ < 1, 

p < \t\^~ < p~ x , with either p* = p,q* — q or p* = q. For n > 1 the crossing 
parameter ry looks like an arbitrary free variable, not related to other parameters 
of the system, but, in fact, it is essentially equivalent to the coupling constant t for 
■Uj-spins. If t — 1, relation (|5.2[) reduces to n-th power of the standard STR. 
Define the i?C„-root system generalization of the F-function: 

(5.7) I(t 1 ,...,t 8 ;t;p,q) = J| T(tjt k ;p,q,t) 

l<j<fe<8 

T(tzfz^;p,q) " ULiWk*? 1 ;P,q) dz 3 



X Kr, 



TT _ *fc ifiVJ TT 

i<M<n M 1 *?™*) M r(^;p 1? ) izj 



where parameters t,ti,...,t$ £ C, satisfy |i|,|tj| < 1, and £ 2n ~ 2 n,_i tj = P 2 Q 2 
constraints. As shown by Rains }R3j . this function obeys the same W(Er) Weyl 
group of symmetries as in the n — 1 case. The key transformation has the form 

(5.8) I(ti, . . .,t s ;t;p,q) = 7(s 1; . . . ,s 8 ;t;p,q), 

where 



s 



,-p-Hj, j = 1,2,3,4 . _ lht 2 t 3 U _ jpqt 



1*1,1^1,^-1 <1. 



s j=P t v J = 5,6,7,8 ' r \ pqt x - n \j t 5 t6t7ts' 

Introduce variables aij by relation t( n-1 " 4 tj — (j>q) 1 ^e 2nlx \ so that the balancing 
condition becomes x j = Then (|5.8p describes the invariance of the integral 

/ with respect to the Weyl reflection 

2(x, v) 

x — » S v {x) — X ; — '—r-v, X, v £ R 8 , 

(v,v) 

where (x, v) = J^t—i x^Vk is the scalar product and the vector v has components 
Vk = 1/2, k = 1,2,3,4, and Vk = —1/2, k — 5,6,7,8. Together with the group 
Ss permuting the parameters Xj, this transformation generates full exceptional 
reflection group W(E?). 

Equality (|5.8[) can be rewritten in the star-star relation form 

f S{u; t, p, q)W( V - a- x, u)W(ri - y, u)W( T , w, u)W(5; z, u) [du] 

J[0,2tt]™ 

(5-9) = R -51——r \ S{u;t,p,q) 

Pt{a + P;w,z) ^[0,2^]" 

x W(/3; x, u)W(a; y, u)VF(?7 — (5; to, u)W(t] — 7; z, u)[du], 
where a + /3 = 7 + 5 and 

In complete analogy with n = 1 case (|1.16[) , (| 1 . 1 7[) , (11.181) , one can obtain two other 
differently looking star-star relations for n > 1 by an iterative application of this 
formula after permutations of parameters. 

Using the matrix integral representations for elliptic hypergeometric integrals, 
111 [SVT] relation (JITS) was shown to describe a new electric-magnetic duality 
between two four-dimensional J\f — 1 supersymmetric Yang-Mills theories with 
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the gauge group G = SP{2n). Namely, the electric theory has the flavor group 
SU(8) x U(l); it contains the vector superfield in the adjoint representation of G, 
one chiral scalar multiplet in the fundamental representations of G and SU(8), and 
the field described by the antisymmetric tensor of the second rank of G. The mag- 
netic theory has the flavor group SU(4:)i x SU(A) r x U(1)b x U(l) and a similar 
set of quantum fields, as well as 2n additional gauge invariant mesonic fields - 
the antisymmetric tensors of S'[/(4)-flavor subgroups. The elliptic Selberg integral 
(|5.1[) describes the confinement phenomenon in the SP(2n) super- Yang-Mills gauge 
theory with 6 chiral supcrficlds in the fundamental and 1 chiral superfield in the 
antisymmetric representations of SP{2n) 1 respectively, — its dual magnetic phase 
contains only a peculiar set of mesonic fields without local gauge symmetry. 

From the present paper point of view relations ()5.2[) and (|5.9j) should have an 
appropriate physical interpretation in the context of discrete integrable models for 
n > 1 similar to the n = 1 case. We described already one possible honeycomb 
lattice model that can be associated with the elliptic Selberg integral. The system 
lying behind relation (|5.9[) resembles the checkerboard Ising model with the contin- 
uous spins. As an elementary cell one has a square with four white vertices (with 
the single component spins x, y, . . . sitting in them) and one black vertex in the 
center (with the n-component spin u sitting in it and the integral taken over its 
values), see figure 4. There are again three differently looking star-star relations for 
n > 1 obtained by repeated application of the same formula (|5.9p in conjugation 
with permutation of parameters, quite similar to the n — 1 case. Equality ()5.2[) 
can be considered then as their reduction to STR. 

We did not discuss in this paper an important physical question of the existence 
of phase transitions in the described models and the spectrum of scaling exponents. 
For clarifying this point it is necessary to single out the temperature like variable 
associated with one of the parameters p or q |BSj and investigate the behavior of the 
partition functions per edge (defined by m(a)'s) when the temperature varies from 
large to small values. Since many known systems with nontrivial phase transitions 
are represented by the limiting cases of the elliptic beta integral STR solutions, there 
are nontrivial critical phenomena. However, their classification requires separate 
analysis and lies beyond the scope of the present work. 

In [SVlj a large number of proven and conjectural evaluation formulas for 
elliptic beta integrals on root systems and their nontrivial symmetry transformation 
analogues for higher order integrals has been listed. Actually, it was conjectured 
that there exist infinitely many such integrals, and for each of them one can expect 
suitable application in the context of solvable models of statistical mechanics and 
other types of integrable systems. 
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Appendix A. The modified (j-gamma function 

The function ^f(u;oj\^ui2) (|3.3p or its various transformed versions are referred 
to in different papers as the double sine function |KLS|, IPTj . the non-compact 
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quantum dilogarithm [B IFKVj fPTl El [BS], the hyperbolic gamma function 
[Rul lB]. or the modified a-gamma function IS9I. 

The functional equations satisfied by j(u;ui) have the form 

{A.L) - - — 1 e "2, ■; r — 1 — e 1 . 

Using a modular transformation for thcta functions one can derive another repre- 
sentation for 7(u;wi,cj2) complementary to Q3.3p : 

The non-compact quantum dilogarithm [F] in the notation of [BMSJ (in [ FKVj 
it was denoted as &b{z)) has the form 

(\ f er 2izw \ dw 

f{z) = exp 



= exp 

where 



4 ,/R+io smh(iufo) sinh(u>/&) / w 

dw 



R+ 



i0 (e wb - \){e w l h -1)J w ' 



u = i(6 + 6 x ) -iz. 
For q — e 27rlfc2 , g = e ^ 27ri / fc2 an d l m (5 2 ) > 0, one can write 

(e 2 ^;^ _ (- 9 1 /2 e 2 7r6 ,. (?)oo 



Therefore, 

V{z)=l{\{b + b- 1 )-iz;b,b- l y 1 - 
The S'b(u) function used in |PTj coincides with 

for u>i = b and u>2 = b^ 1 , and another similar function of |PT| is 
G b {u) = e-T5( 3+b2 +^ 2 ) 7 (u;6,6- 1 ). 

The 7(z)-function used in [Vi coincides with the infinite products ratio on the 
right-hand side of (|A.2[) for z = u>2U and r = ui%/u)2- 

For Re(wi), Re(oj2) > and < Re(w) < Re(wi+W2) the function 7^ 2 )(m; oji, LU2) 
has the following integral representation 

(A.3) 7 (2) (u;wi,w 2 ) = exp ( -PV 



e ux dx 



(1 - e" lX )(l - e^ x ) x 

where the principal value of the integral means PV J R = 2~ 1 (f R+ . Q + J R _ i0 )- Using 

the fact that PV J R dx/x k = for k > 1, one can write 

(r) ( f°° t sinh(2it — u>i — ui2)x 2u — uj± — u>2\\ dx 

7^(u;wi,w 2 ) = exp - / I r-— ; J — . 

\ J V 2 smn(cjix) smn(w 2 ^J ZbJiuj 2 x Ijx 

Comparing this expression with the hyperbolic gamma function Gh{z;uj) defined 
in [Ruj . one can see that 

G h (z;u) = 7 ( 2 )(i( Wl +uj 2 ) -iz;uj, Re(w 1 ),Re(w 2 ) > 0. 
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Changing in (|A.3|) the sign of the integration variable x — > — x and simultaneously 

ujk — > —LOk, u — > — u, we find that 



(A.4) 7 < 2 > (u; ui , wa) = exp ( +PV 



dx 



(f - e^ x )(l - e^ x ) x J ' 

where Re(wi), Re(w2) < and Re(wi + UI2) < Re(u) < 0. 

The double sine function is dehned as 5*2 (u; w) = l/ 7 ( 2 )(it; cj) and its properties 
were described in detail in the Appendix of [KLS] ■ For the 7^ 2 )-function we have 

(2)f \* (2)f * * *\ (2)/ W l+ w 2 , \ , 

7 V ; (u;wi,w 2 ) = T >{u ; u 1 , oj 2 ), V ( ^ ± tt; Wi, w 2 ) = 1, 

and 7^ 2 )(au; awi, 01^2) = 7 < - 2 '(u; wi,W2) for arbitrary complex o ^ 0. After such a 
rescaling in (IA.4I) with a = 2iri one gets the definition of the hyperbolic gamma 
function given in |R2j . 

The asymptotics we are interested in for Im(o>i / 102) > have the form 

lim eT Bl ' , (" ;u )'y( 2 )(tj; w) = 1, for argwi < argit < argo>2 + n, 

u—toc 

lim e _ ~ B2 ' 2< - u; '* ; ' l 7 < ' 2 - ) (u; u) = 1, for argwi — 7r < argu < argW2- 

Appendix B. General multiple gamma functions 

Barnes introduced a multiple zeta function as the following m-fold series [Bar 
00 ^ 

(m(s,u;uj) = ^2 i , ma . = mcJi + . . . + n m w m , 

nx,...,n m =0 + 

where s, u G C. This series converges for Re(s) > m under the condition that all u>j 
lie in one half-plane dehned by a line passing through zero. Because of the latter 
requirement, the sequences n\U)\ + . . . + n m uj m do not have accumulation points on 
the hnite plane for any rij —¥ +00. It is convenient to assume for definiteness that 
Re(wj) > 0. 

The function £ m (s,w;w) satishes equations 

(B.l) C, m {s,u + u)j\ui) - £ m {s,u;u)) = -C m -i(s,tt;u;(j)), j = l,...,m, 

where ui(j) = (oj\, . . . ,ojj-x,ojj±x, . . . ,w m ) and £o(s,u;u>) = u~ s . The Barnes mul- 
tiple gamma function is defined by the equality 



s=0 



T m (u;u>) = cxp(d( m (s, u;uj)/ds) 
It satisfies finite difference equations 

(B.2) T m {u-\-u)j\ui) = t~ — 7-Trr m (w;w), j = l,...,m, 

where Tq(u;ui) := u" 1 . 

The multiple sine-function is defined as 

S m (u;w) = — ^ k=1 ' ' 



T m (u;u) 

It is more convenient to work with the hyperbolic gamma function 

1 i - m \u-Lu) = S m {u-Lu)<-- 1 ^ 1 
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satisfying the equations 

7 M (u + Uj ; U ) = 7 (m - 1} (u; w(j)) 7 (m) (u; w), j = 1, . . . , m. 

Note that the elliptic gamma function can be written as a special combination of 
four Barnes gamma functions of the third order |S9j , and similar relations are valid 
for higher order elliptic gamma functions used in the present paper. 
One can derive the integral representation [N] 

(m) - - ' — r eUX dx 



7 lmJ (u;o;) = exp -PV 



7T1 

exp ( rB m , m («;w) 



nr=i(^^-i) x 

e ux dx 



io nr=i(^-i)' 



i! 

-io nr=i(^-i)' 

where Re(ojfe) > and < Re(tt) < Re(X)fcLi w fe) an d B m ^ m are multiple Bernoulli 



exp ( —B mm (u;w) - \ 

id < Re(tt) < Re(£™ =1 w 
polynomials defined by the generating function 



n 



(b.3) nr=i( e " fc!c - 1) = o gm - n(u; • • • >fa,m) ^r- 

Infinite product representations for these functions have been derived in [N] , 
In particular, for [p|, |g| < 1 and |r| > 1 we have 

r 3 ), , .-**,.,(,.;«) fj (l-e 2 ^Mg 3+ i r -( fc+ i))(i- e 2^ Vgfc ) 

' \ ' J 11 ^ _ e 2mu/u s fjj+l r -k 

j,k=0 ' ^ 

which is used in the main text after the reduction uj^ = 2{oj\ + CJ2 ) (or p = q 2 , 

r — p _ g-7riw 2 /(wi+W2)^ 

The functions m(a) (|4.3p . (|4.6[) . and (|5.6[) defining the free energy per edge as 
described in the main part of the paper are related to particular cases of the Lerch 
type generalization of the Barnes zeta-function: 

00 n m p nh 

( m (s,u;f3;uj) = (m + 0? ' ^ = n^i + . . . + n m uj m , 

ni^....nm—0 

converging for all \Pk\ < 1, or Re(s) > m and \Pk\ — 1 (provided the same con- 
straints on ujj are valid as in the plain Barnes case). The univariate case, i.e. the 
proper Lerch zeta-function, is described, e.g., in |WWj . 

The function ( m (s,u; f3;ui) satisfies the following set of finite difference equa- 
tions 
(B.4) 

pjC m (s,u + Wj;P;cj) - £ m (s, u; p; u) = -Cm-i(s, u; P(j); j = l,...,m, 

where = . . . , Uj-i, u j+ i, . . . , w m ), P(j) = (Pi, • - • , Pj-i, Pj+i, ■ • - , Pm), 

and Co(s, u; P; w) = u~ s . 

Similar to the Barnes case, one can easily derive the integral representations 

1 f°° t s ~ 1 e~ ut 
( m (s,u;P;uj) = ——- m — -rr dt 



,rri - a) f {-ty-'e-^ 
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where Ch is the Hankel contour encircling the half-line [0,oo) counterclockwise, 
and using them analytically continue £ m -function in s and /3fe to different regions 
of parameters. The /^-deformation of the Barnes multiple gamma function defined 
as r m (u; /3; ui) — exp(dCm (s,u;/3;w)/ds)\ s=Q satisfies the finite difference equations 

(B.5) T m (u + ojjifcojfj = - T m (u;j3;u), j = l,...,m, 

where T (u; ft; oj) := it -1 . 

When ftk are primitive roots of unity, f3^ k = 1, rife = 2, 3, . . . , it is possible to 
rewrite £ m (s,it; /3; ui) as linear combinations of the standard Barnes zeta functions. 
It follows from the simple identity 

1 = n i =0,2,...,n t -l(l-flU) 

1 - /3 k z 1 - z n " 

This allows expressing the functions like (|4.5p as linear combinations of the standard 
Barnes gamma functions, which was used in the construction of infinite product 
representations of the functions m(a) (|4.3p . (|4.6I) . and (I5.6[) . In particular, function 
(|4.5|) is emerging from the m = 3 case with the choice /3i = f3 2 = 1, /?3 = — 1. 
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